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Abstract 

In  this  paper,  we  use  a generic  and  general  variational  method  to  obtain 
solutions  to  the  flow  of  generalized  Newtonian  fluids  through  circular  pipes  and 
plane  slits.  The  new  method  is  not  based  on  the  use  of  the  Euler-Lagrange 
variational  principle  and  hence  it  is  totally  independent  of  our  previous 
approach  which  is  based  on  this  principle.  Instead,  the  method  applies 
a very  generic  and  general  optimization  approach  which  can  be  justified 
by  the  Dirichlet  principle  although  this  is  not  the  only  possible  theoretical 
justification.  The  results  that  were  obtained  from  the  new  method  using 
nine  types  of  fluid  are  in  total  agreement,  within  certain  restrictions,  with 
the  results  obtained  from  the  traditional  methods  of  fluid  mechanics  as  well 
as  the  results  obtained  from  the  previous  variational  approach.  In  addition 
to  being  a useful  method  in  its  own  for  resolving  the  flow  field  in  circular 
pipes  and  plane  slits,  the  new  variational  method  lends  more  support  to 
the  old  variational  method  as  well  as  for  the  use  of  variational  principles  in 
general  to  resolve  the  flow  of  generalized  Newtonian  fluids  and  obtain  all  the 
quantities  of  the  flow  field  which  include  shear  stress,  local  viscosity,  rate  of 
strain,  speed  profile  and  volumetric  flow  rate.  The  theoretical  basis  of  the 
new  variational  method,  which  rests  on  the  use  of  the  Dirichlet  principle,  also 
provides  theoretical  support  to  the  former  variational  method. 

Keywords:  fluid  dynamics;  rheology;  variational  method;  Dirichlet  principle; 
pipe  flow;  slit  flow;  Newtonian;  power  law;  Ellis;  Ree-Eyring;  Carreau;  Cross; 
Bingham;  Herschel-Bulkley;  Casson. 
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1 Introduction 


The  flow  through  circular  pipes  and  plane  slits  has  many  applications  in  physical 
and  biological  sciences  and  engineering  and  hence  it  has  been  investigated  in  the  past 
by  many  researchers  (e.g.  [1-12])  using  various  methods  of  fluid  dynamics.  Recently 
we  proposed  the  use  of  Euler-Lagrange  variational  principle  [13]  to  resolve  the  flow 
of  generalized  Newtonian  fluids  through  circular  pipes.  The  method  is  based  on 
minimizing  the  total  stress  in  the  flow  conduit  in  the  sense  of  minimizing  the  stress 
profile  in  the  velocity- varying  dimension.  This  attempt  was  later  expanded  and 
supported  by  other  investigations  [14-16]  where  the  method  was  successfully  applied 
to  more  types  of  fluid  and  another  type  of  geometry,  namely  the  plane  slit  conduit. 

Despite  the  success  of  this  method  in  describing  the  flow  of  several  fluid  models 
and  conduit  types,  it  has  not  been  proven  in  general  by  a formal  mathematical 
argument  that  justifies  the  universal  applicability  of  the  variational  method  and 
the  principle  on  which  it  relies.  Certain  mathematical  technicalities  may  also  be 
disputed  and  hence  it  is  desirable  to  fortify  the  method  by  a more  generic  and 
general  variational  approach  that  is  more  safe  from  such  disputes  and  pitfalls.  The 
present  investigation  tries  to  do  so  by  using  a very  basic  and  general  variational 
approach  where  the  flow  field  in  the  conduit  is  resolved  through  the  application  of 
a generic  optimization  technique  to  a stress  functional.  The  theoretical  justification 
of  this  functional  can  be  obtained  from  the  Dirichlet  principle  although  it  can  also 
be  justified  by  other  theoretical  foundations  based  on  purely  physical  arguments. 

The  plan  for  this  paper  is  that  in  section  § 2 we  present  a general  description  of 
the  proposed  method  and  its  theoretical  background.  This  is  followed  in  section  § 
3 by  discussing  practical  issues  about  the  implementation  of  this  method  and  the 
presentation  of  sample  results  that  were  obtained  from  this  method  with  comparison 
to  similar  results  obtained  from  the  previous  methods  which  include  the  classical 
methods  of  fluid  mechanics  and  the  former  variational  approach.  The  paper  is 
finalized  in  section  § 4 with  general  discussions  and  summarization  of  the  main 
achievements  of  the  present  study.  As  a matter  of  convenience,  we  label  the  former 
variational  method  which  is  based  on  the  Euler-Lagrange  principle  with  EL  and 
the  new  variational  method  which  is  based  on  the  Dirichlet  principle  with  DM. 
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2 Method 


We  first  state  our  assumptions  about  the  flow,  fluid  and  conduit  which  are  adopted 
in  the  present  study.  We  assume  a laminar,  isothermal,  incompressible,  steady, 
pressure-driven,  fully-developed  flow  of  a time-independent,  purely-viscous  fluid 
that  can  be  described  by  the  generalized  Newtonian  fluid  model,  that  is 


where  the  viscosity,  fi,  and  stress,  r,  depend  only  on  the  contemporary  rate  of  strain, 
7,  and  hence  the  fluid  has  no  memory  of  its  deformation  past.  In  this  formulation 
we  ignore  all  non  deformation-related  dependencies  of  the  viscosity  and  stress  due 
to  other  physical  factors  like  temperature  and  pressure.  In  fact  we  consider  only  the 
shearing  effects  since  the  effects  of  other  forms  of  deformation,  such  as  extensional, 
are  presumed  insignificant  which  is  well  justified  for  the  presumed  state  of  flow, 
fluid  and  conduit.  Edge  effects  at  the  entry  and  exit  of  the  conduit,  as  well  as 
external  body  forces,  are  also  regarded  negligible. 

Concerning  the  conduit,  we  use  circular  pipe  and  plane  slit  geometries,  which 
are  depicted  in  Figures  1 and  2,  where  the  pipe  is  assumed  straight  with  a cross 
section  that  is  uniform  in  shape  and  size  while  the  slit  is  assumed  straight  long 
and  thin  with  a uniform  cross  section.  In  both  cases  we  assume  rigid  mechanical 
characteristics  of  the  conduit  wall  as  opposite  to  being  deformable  such  as  having 
elastic  or  viscoelastic  mechanical  properties.  It  is  also  assumed  that  the  slit  is 
positioned  symmetrically  in  its  thickness  dimension,  z,  with  respect  to  the  plane 
£ = 0 as  depicted  in  Figure  2. 

Regarding  the  boundary  conditions,  we  assume  no-slip  at  the  conduit  wall,  where 
the  fluid  interfaces  the  solid  [17],  with  the  flow  speed  profile  having  a stationary 
derivative  point  at  the  symmetry  center  line  of  the  pipe  and  symmetry  center  plane 
of  the  slit  which  means  zero  stress  and  rate  of  strain  at  these  loci.  As  for  the 
viscoplastic  fluids,  this  stationary  region  expands  to  include  all  the  points  at  the 
forefront  of  the  flow  profile  whose  stress  falls  below  the  yield  stress,  as  will  be 
discussed  further  in  the  coming  sections. 

Now,  for  the  generalized  Newtonian  fluids  that  satisfy  the  above  assumptions, 
the  momentum  equation  in  one  dimension  is  reduced  to 


T = VI 


(1) 


ds 


(2) 
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Figure  2:  A schematic  of  the  plane  slit  geometry  used  in  this  investigation. 


where  s is  a spatial  coordinate  that  represents  r for  pipes  and  z for  slits,  and  G is  a 
constant.  If  we  differentiate  the  last  equation  with  respect  to  s we  get 


(3) 


which  is  a one- variable  Laplace  equation  in  one  dimension.  According  to  the  Dirichlet 
principle,  the  solution  of  this  equation  is  a minimize!’  of  the  following  functional 


and  vice  versa 
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(4) 


D = 


ds 


where  D is  the  Dirichlet  functional  and  Q is  a spatial  domain  in  this  formulation. 
On  using  Equation  2,  substituting  and  simplifying,  the  functional  can  be  reduced  to 


D = \G\  j dr  (5) 

where  T is  the  new  domain  in  this  formulation.  The  last  equation  demonstrates  that 
the  solution  of  this  problem  is  a minimizer  (and  vice  versa)  of  the  total  stress  in  the 
sense  that  have  been  given  previously  and  hence  it  establishes  the  foundation  of  our 
former  variational  approach,  EL,  which  is  based  on  the  use  of  the  Euler-Lagrange 
principle.  It  can  also  provide  a theoretical  foundation  for  new  method  or  methods. 

Apart  from  its  theoretical  aspects,  the  Dirichlet  principle  can  be  used  practically 
as  a basis  for  another  variational  method,  DM,  that  can  be  employed  to  obtain 
flow  solutions  and  verify  the  solutions  obtained  by  other  methods  including  the 
EL  method.  The  DM  method  in  practical  terms  is  based  on  finding  the  shear 
stress  solution  in  conduits  by  minimizing  the  above  functional.  This  functional  is 
discretized  and  minimized  numerically  using  an  optimization  algorithm  subject  to 
the  boundary  conditions  at  the  conduit  wall  and  conduit  center.  For  pipes  these 
boundary  conditions  are  respectively 


Ad  — A? 


RAp 
2 L 


while  for  slits  they  are 


and 


An  = 0 


(6) 


tw  = tb=  ——  and  rm  = 0 (7) 

where  rw  is  the  shear  stress  at  the  conduit  wall  which  is  equivalent  to  Tr  for  pipes 
and  to  tb  for  slits,  rm  is  the  stress  at  the  conduit  center  line  or  plane,  R is  the 
pipe  radius,  B is  the  slit  half  thickness,  and  L is  the  conduit  length  across  which  a 
pressure  drop  A p is  exerted. 

The  numerically  obtained  solution,  r(s),  is  then  used  in  conjunction  with  the 
rheological  constitutive  relations,  as  given  in  Table  1 for  the  models  considered  in 
this  study,  applied  to  the  generalized  Newtonian  fluid  equation  to  find  y(s)  either 
explicitly  or  implicitly  through  the  use  of  a simple  numerical  solver  like  a bisection 
solver.  As  for  the  viscoplastic  fluids,  a zero  stress  is  applied  to  all  the  points  at 
the  forefront  of  the  flow  profile  whose  shear  stress  falls  below  the  yield  stress  value 
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during  the  optimization  process.  As  indicated  earlier,  this  is  an  extension  to  the 
second  boundary  condition  at  the  conduit  center  to  include  a plane  region  in  the 
forefront  of  the  speed  profile  and  hence  it  does  not  compromise  the  optimization 
condition  and  the  underlying  variational  principle. 

The  obtained  y(s)  is  then  integrated  numerically  with  respect  to  s to  find  the 
flow  speed,  w(s),  where  the  no-slip  boundary  condition  at  the  conduit  wall  is  used 
to  provide  an  initial  value,  v = 0,  that  is  incremented  on  moving  inward  from  the 
conduit  wall  toward  the  conduit  center  during  the  integration  process.  This  is 
followed  by  integrating  u(s)  numerically  with  respect  to  the  conduit  cross  sectional 
area  normal  to  the  flow  direction  to  obtain  the  volumetric  flow  rate.  During  these 
successive  integration  processes,  the  boundary  conditions  at  the  conduit  wall  and 
center,  which  are  based  on  the  zero  speed  and  zero  stress  respectively,  are  used. 

In  Table  1 the  rheological  constitutive  relations  for  the  nine  fluid  models  which 
are  employed  in  this  study  are  presented,  while  in  Tables  2 and  3 the  analytical 
relations  that  correlate  the  flow  rate,  Q,  to  the  applied  pressure  drop,  A p,  for  the 
flow  in  pipes  and  slits  respectively  are  given.  Most  of  these  expressions  can  be 
found  in  the  classic  literature  of  rheology  and  fluid  dynamics  (e.g.  [18,  19])  while 
the  rest  can  be  obtained  with  their  derivation  from  [14-16].  Regarding  the  Carreau 
and  Cross  fluids,  the  factors,  which  are  included  in  their  Q expressions  and 
represent  definite  integral  expressions,  are  given  by 


'PiCa 


+ 


h3  [3A4  (3n/2  + 5 n'  + 2)  74  - 3n'A27|  + 2]  (1  + A2y 2Rfn'/2 
3A4  (9n'2  + 18n'  + 8) 

Hi52  [A4  (2 n'2  + 5 n'  + 3)  7^  - n'\2^2R  + 1]  (1  + \2r)R)n 
2A4  (n'  + 1)  in'  + 2) 

tfd  [A4  {n'2  + 5 n'  + 6)  7 4R-  n'\2^2R  + 2]  (1  + X2^2R)n'/2  + 


A4  (ra'  + 2)  (n'  + 4) 
253  m52 


+ 


+ 


2^5 


3A4  (9n'2  + 18n'  + 8)  2A4  (n'  + 1)  (V  + 2)  A4  {n>  + 2)  {n>  + 4) 


(8) 


'p,Cr 


{2(53  [— m (2 f2  + 5/  + 3)  + 4 g2  + 2m2]  + 12 mS2^  (m  — g)  + 12 m25n2g2  + 3 m2^i3g3]  jR 

12  m2g3 

{53  (m2  - 6m.  + 8)  + 3md2fM  (m  - 4)  + 3 m2<)>2}  2F\  (l,  m±A- -/) 

12m2  1 j 
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*s,Ca 


n'52^B  [2F1  (§,  1 - n'\  §;  -A27|)  - 2F\  (1,  -ra';  |;  — A2^)] 


A2 


+ 


(1  + wQ  82 7I  2F1  (§,  -w';  f ; -A27|) 

3 

n'&miB  [2-P1  (§, 1 - f;  -A27l)  - 2-Pi  (§,  -y;  |;  -A27b)] 

A2 

(2  + n')  8na%  2^1  (|,  - y;  -A27l)  + All 


(10) 


and 


4,Cr  — 


[352  (■ m-g ) - {52  (m  - 3)  + 2m<5/^}  #2  2Fi  (l,  A;  1 + A;  -/)  + 6 mSfng  + 2 mrfg2]  7 1 


6 m g 2 


(11) 


where,  in  these  expressions, 


5 = {no-  in) , ri  = {n-  1) , / = Am7™,  g = 1 + /,  (12) 

and  2 -Pi  is  the  hypergeometric  function  of  the  given  arguments  with  its  real  part 
being  used  in  the  evaluation  of  / factors.  Moreover 


HrIr  = i~R  and  iibIb  = tb  (13) 

where  Tr  and  Tb  are  given  by  Equations  6 and  7 respectively,  with  7 r and  7 b being 
obtained  numerically  from  the  above  implicit  relations,  as  explained  in  [16]. 
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Table  1:  The  constitutive  relations  for  the  nine  fluid  models  used  in  this  investigation. 
The  meaning  of  the  symbols  are  given  in  Nomenclature  § 5. 


Model 

Constitutive  Relation 

Newtonian 

r = 

ho7 

Power  Law 

r = 

/ \ a— 1 

Ellis 

h = 

he 

i 

+ 

Ree-Eyring 

T = 

rcarcsinh 

Carreau 

h = 

hi  + (ho  — hi)  (1 

+ A272)(”-1)/2 

Cross 

h = 

II  1 A^o  ~^i 

f'l  T 

Bingham 

T = 

C"7  + To 

Herschel-Bulkley 

T = 

Cl 

n + r0 

Casson 

T1/ 2 

= (A'7)1/2  + rP 

Table  2:  The  volumetric  flow  rate,  Q , for  the  pipe  flow  of  the  nine  fluid  models 
used  in  this  investigation.  The  symbols  are  given  in  Nomenclature  § 5. 
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Table  3:  The  volumetric  flow  rate,  Q,  for  the  slit  flow  of  the  nine  fluid  models  used 
in  this  investigation.  The  symbols  are  given  in  Nomenclature  § 5. 


Model 
Newtonian 
Power  Law 

Ellis 

Ree-Eyring 

Carreau 

Cross 

Bingham 

Herschel-Bulkley 

Casson 


Q 


2WB3Ap 


Sfi0L 

2WB2n 

2n+l 

V 

2WB2 

Me 

3 

2tVr3 

' B 

Mr 

\TB 

2 WB2Is,c* 

2WB2IStCl 

BAp 

kL 


hr1 


(o+2)r; 

tb  cosh  ( ^ ) - tc  sinli  ( ^ 


2W  1 

B_ 

C'  \ 

2 W 

( — 

Vc 

VTB 

2 W 1 

'b_' 

~K  \ 

,TB  j 

2 r 


A 

3 


To  Tl 


+ A 
T 6 


2 r 


n(nTQ+nTB+TB)(TB-To)1+1/n 


r3 
I B_ 

3 


(2n2+3n+l) 

4P^4/2  i _ zl 

5 “r  2 30 
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3 Implementation  and  Results 


The  optimization  method,  as  described  in  the  last  section,  was  implemented  in  a 
computer  code  using  five  numerical  optimization  algorithms:  three  deterministic 
which  are  Conjugate  Gradient,  Quasi-Newton  and  Nelder-Mead  [20];  and  two 
stochastic  which  are  the  Stochastic  Global  algorithm  of  Boender  et  al  [21]1  and  a 
generic  simulated  annealing  algorithm  [22],  These  five  algorithms  produce  similar 
solutions  with  different  levels  of  accuracy  and  convergence  rate.  The  sample 
results  presented  in  this  paper  are  obtained  mostly  from  the  Stochastic  Global 
algorithm  which  is  overwhelmingly  the  most  accurate  and  reliable  one.  In  addition, 
standard  numerical  integration  and  bisection  solution  techniques,  as  well  as  standard 
algorithms  for  evaluating  complicated  functions  like  hypergeometric,  were  employed. 

The  newly  proposed  variational  method  was  then  employed  to  obtain  solutions 
for  the  flow  of  nine  types  of  fluid  through  pipes  and  slits.  The  nine  types  of  fluid 
are:  Newtonian,  power  law,  Ellis,  Ree-Eyring,  Carreau,  Cross,  Bingham,  Herschel- 
Bulkley  and  Casson.  The  results  obtained  from  the  new  method  using  wide  ranges 
of  fluid  and  conduit  parameters  were  thoroughly  compared  to  the  results  obtained 
from  the  traditional  methods  of  fluid  mechanics  and  the  former  variational  method. 

In  all  the  investigated  cases  the  three  methods  produced  very  similar  results 
within  acceptable  error  margins.  One  exception  is  the  Ellis  model  for  which  EL  has 
not  been  formulated  and  implemented.  Another  exception  is  the  viscoplastic  fluids 
where  the  EL  method  differs  significantly  when  the  yield  stress  value  is  high.  This 
is  justified  by  the  fact  that  the  EL  method  was  formulated  and  implemented  for  the 
non- viscoplastic  fluids  specifically  and  hence,  as  we  demonstrated  in  our  previous 
investigations  [13-15],  it  is  just  an  approximation  for  the  viscoplastic  fluids  which  is 
a good  one  only  when  the  yield  stress  value  is  low.  However,  we  believe  that  even 
the  EL  method  can  be  re-formulated  and  re-implemented  to  include  viscoplastic 
fluids,  and  hence  it  can  produce  similar  results  to  the  other  methods  even  for  the 
high  yield  stress  fluids,  although  we  did  not  make  any  effort  in  the  current  study  to 
do  so. 

A representative  sample  of  these  results  obtained  from  the  three  methods  are 
presented  in  Figures  3,  4,  5 and  6 where  the  fluid  and  conduit  parameters  of  these 
examples  are  given  in  Tables  4 and  5.  In  these  figures,  the  analytical  solution  is 
represented  by  the  solid  line  while  the  solutions  from  the  two  variational  methods 
are  represented  by  the  circles.  The  exception,  as  indicated  earlier,  is  the  Ellis  and 

^ee  also:  http://jblevins.org/mirror/amiller/global.txt  webpage. 
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viscoplastic  fluids  where  the  circles  represent  only  the  DM  method.  The  reason  for 
combining  the  solutions  of  the  two  variational  methods  in  a single  representation  is 
that  they  produce  very  similar  results  and  hence  there  is  no  point  in  plotting  them 
separately. 

As  seen  in  these  examples,  the  two  variational  solutions  agree  very  well  with  the 
analytical  solutions.  The  minor  departure  in  some  cases  between  the  two  variational 
methods  on  one  hand  and  the  analytical  on  the  other  is  due  mainly  to  the  nature 
of  the  variational  methods  as  they  heavily  rely  on  numerical  techniques,  mainly 
bisection  solvers  and  numerical  integration,  which  is  not  the  case  with  the  analytical 
solutions  since  they  are  evaluated  directly.  As  indicated  above,  unlike  the  EL 
approach  which  is  a good  approximation  for  the  viscoplastic  fluids  only  if  their 
yield  stress  is  low,  the  DM  approach  produces  ‘exact’  solutions,  considering  the 
numerical  errors  introduced  by  the  heavy  use  of  numerical  techniques,  even  for  the 
fluids  with  high  yield  stress. 
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(a)  Newtonian  (b)  Power  Law  (c)  Ellis 


(d)  Ree-Eyring  (e)  Carreau  (f)  Cross 


(g)  Bingham  (h)  Herschel-Bulkley  (i)  Casson 


Figure  3:  Comparing  the  analytical  solution  (solid  line)  to  the  variational  solutions 
(circles)  of  7 in  s-1  (vertical  axis)  versus  r in  m (horizontal  axis)  for  the  flow  of  the 
nine  fluid  models  in  pipes.  The  EL  solutions  are  not  represented  for  the  Ellis  and 
viscoplastic  fluids.  The  pipe  and  fluid  parameters  are  given  in  Table  4,  where  in  all 
cases  A p = 500  Pa. 
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(a)  Newtonian  (b)  Power  Law  (c)  Ellis 


(d)  Ree-Eyring  (e)  Carreau  (f)  Cross 


Figure  4:  Comparing  the  analytical  solution  (solid  line)  to  the  variational  solutions 
(circles)  of  7 in  s-1  (vertical  axis)  versus  z in  m (horizontal  axis)  for  the  flow  of  the 
nine  fluid  models  in  slits.  The  EL  solutions  are  not  represented  for  the  Ellis  and 
viscoplastic  fluids.  The  slit  and  fluid  parameters  are  given  in  Table  5,  where  in  all 
cases  A p = 700  Pa. 
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(a)  Newtonian  (b)  Power  Law  (c)  Ellis 


Figure  5:  Comparing  the  analytical  solution  (solid  line)  to  the  variational  solutions 
(circles)  of  Q in  m3.s-1  (vertical  axis)  versus  A p in  Pa  (horizontal  axis)  for  the  flow 
of  the  nine  fluid  models  in  pipes.  The  EL  solutions  are  not  represented  for  the  Ellis 
and  viscoplastic  fluids.  The  pipe  and  fluid  parameters  are  given  in  Table  4. 
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(a)  Newtonian  (b)  Power  Law  (c)  Ellis 


(d)  Ree-Eyring  (e)  Carreau  (f)  Cross 


Figure  6:  Comparing  the  analytical  solution  (solid  line)  to  the  variational  solutions 
(circles)  of  Q in  m3.s-1  (vertical  axis)  versus  A p in  Pa  (horizontal  axis)  for  the  flow 
of  the  nine  fluid  models  in  slits.  The  EL  solutions  are  not  represented  for  the  Ellis 
and  viscoplastic  fluids.  The  slit  and  fluid  parameters  are  given  in  Table  5. 
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Table  4:  Fluid  and  pipe  parameters  for  the  examples  of  Figures  3 and  5.  SI  units 
apply  to  all  dimensional  quantities  as  given  in  Nomenclature  § 5. 


Model 

Fluid  Properties 

R 

L 

Newtonian 

Ho  = 0.025 

0.08 

0.65 

Power  Law 

A;  = 0.033,  n — 1.15 

0.01 

0.095 

Ellis 

He  = 1.42,  Th  = 15,  a = 3.3 

0.04 

0.15 

Ree-Eyring 

Hr  = 0.02,  tc  = 300 

0.03 

0.8 

Carreau 

/i0  = 0.1,  Hi  = 0.008,  A = 1.2,  n = 0.65 

0.03 

0.45 

Cross 

Ho  = 0.15,  Hi  = 0.005,  A = 7.9,  m = 0.8 

0.012 

0.15 

Bingham 

C'  = 0.037,  To  = 2.5 

0.01 

0.1 

Herschel-Bulkley 

C = 0.47,  r0  = 5.0,  n = 0.75 

0.04 

0.7 

Casson 

K = 0.75,  r0  = 3.0 

0.09 

1.33 

Table  5:  Fluid  and  slit  parameters  for  the  examples  of  Figures  4 and  6,  where 
in  all  cases  W = 1.0  m.  SI  units  apply  to  all  dimensional  quantities  as  given  in 
Nomenclature  § 5. 


Model 

Fluid  Properties 

B 

L 

Newtonian 

Ho  = 0.13 

0.02 

1.2 

Power  Law 

k = 0.073,  n = 0.67 

0.01 

0.95 

Ellis 

He  = 0.049,  Th  = 5.0,  a = 2.9 

0.011 

1.95 

Ree-Eyring 

Hr  = 0.57,  tc  = 75 

0.05 

11.5 

Carreau 

Ho  = 0.32,  Hi  = 0.096,  A = 0.75,  n = 0.85 

0.023 

0.75 

Cross 

Ho  = 0.015,  Hi  = 0.007,  A = 2.56,  m = 0.73 

0.035 

2.13 

Bingham 

C'  = 0.48,  r0  = 4.3 

0.014 

1.03 

Herschel-Bulkley 

C = 0.03,  r0  = 5.2,  n = 1.45 

0.035 

3.09 

Casson 

K = 0.12,  r0  = 2.15 

0.017 

2.33 
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4 Conclusions 


In  this  paper  we  presented  a variational  approach  for  finding  the  flow  solutions 
in  one  dimensional  flow  that  applies  easily  to  circular  pipes  and  plane  slits.  The 
method,  which  is  demonstrated  using  nine  types  of  fluid,  can  be  employed  to  obtain 
all  the  required  flow  parameters  which  include  shear  stress,  local  viscosity,  shear  rate, 
speed  profile  and  volumetric  flow  rate.  We  also  presented,  through  the  application 
of  the  Dirichlet  principle,  a theoretical  justification  for  the  application  of  minimizing 
the  total  stress  profile  as  a basis  for  our  variational  approaches  including  the  EL 
method. 

Thorough  comparisons  were  made  both  to  the  analytical  solutions  obtained  from 
the  traditional  methods  of  fluid  dynamics  and  to  the  analytical  or  semi  analytical 
solutions  obtained  from  the  EL  method.  In  all  cases,  the  three  methods  produced 
very  close  results  where  the  differences  can  be  explained  by  the  numerical  errors 
introduced  by  heavy  use  of  numerical  methods  like  numerical  integration,  bisection 
solvers  and  numeric  evaluation  of  complicated  functions.  The  exception  is  the 
viscoplastic  fluids  for  which  the  EL  method  cannot  provide  reliable  solutions  when 
the  yield  stress  is  high  due  to  the  particular  formulation  and  implementation  of 
this  method  which  is  restricted  to  non-viscoplastic  fluids.  The  EL  method  also  has 
not  been  formulated  and  implemented  in  the  past  for  the  Ellis  fluid  and  we  did  not 
make  any  effort  to  do  so  in  the  current  study. 

Apart  from  being  useful  on  its  own  for  resolving  the  flow  field  in  the  given 
conduits  and  obtaining  all  the  required  parameters,  the  proposed  DM  method 
adds  more  support  to  the  previous  EL  approach,  which  is  based  on  applying  the 
Euler-Lagrange  variational  principle,  as  it  confirms  the  results  obtained  from  the  EL 
method  and  provides  a theoretical  foundation  for  it.  Although  the  new  method  may 
not  be  conceptually  identical  to  the  previous  one,  it  should  still  lend  support  to  the 
previous  one  not  only  because  the  two  differently-formulated  variational  methods 
produce  similar  results  but  also  because  they  are  both  based  on  similar  variational 
principles.  From  a procedural  point  of  view,  the  two  methods  are  equivalent  because 
what  is  done  in  DM  numerically  is  done  in  EL,  as  presented  in  our  previous  studies, 
either  analytically  or  partly  analytically  and  partly  numerically.  It  has  also  been 
shown  that  the  theoretical  foundation  of  the  DM  method,  represented  by  the 
Dirichlet  principle,  endorses  the  particular  formulation  of  stress  minimization  which 
EL  rests  upon. 

There  is  an  obvious  theoretical  value  of  the  new  DM  variational  method  which 
is  more  important  than  its  practical  value  which  may  be  insignificant  for  the 
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investigated  fluid  types  and  conduit  geometries  due  to  the  availability  of  the 
presented  analytical  and  numerical  flow  solutions  from  other  methods.  If  the 
proposed  variational  principle  enjoys  a general  applicability,  which  is  yet  to  be 
established  beyond  the  one  dimensional  flow  through  pipes  and  slits,  the  method 
may  also  have  a significant  practical  value  for  the  flow  systems  which  are  more 
complicated  than  the  flow  systems  in  pipes  and  slits  where  the  DM  variational 
method  may  provide  solutions  that  other  methods  might  fail  to  provide,  or  the  DM 
method  may  require  less  effort  to  obtain  the  solutions  than  the  effort  required  by 
the  other  methods. 
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5 Nomenclature 

a indicia!  parameter  in  Ellis  model 
7 rate  of  shear  strain  (s_1) 

7 b rate  of  shear  strain  at  slit  wall  (s-1) 

7ij  rate  of  shear  strain  at  pipe  wall  (s_1) 
jw  rate  of  shear  strain  at  conduit  wall  (s-1) 

<5  ho  - hi  (Pa.s) 

A characteristic  time  constant  in  Carreau  and  Cross  models  (s) 

/i  fluid  shear  viscosity  (Pa.s) 

ho  zero-shear  viscosity  in  Carreau  and  Cross  models  (Pa.s) 
fxe  low-shear  viscosity  in  Ellis  model  (Pa.s) 

Hi  infinite-shear  viscosity  in  Carreau  and  Cross  models  (Pa.s) 
ho  Newtonian  viscosity  (Pa.s) 

hr  characteristic  viscosity  in  Ree-Eyring  model  (Pa.s) 
r shear  stress  (Pa) 

To  yield  stress  in  Bingham,  Herschel-Bulkley  and  Casson  models  (Pa) 
Tb  shear  stress  at  slit  wall  (Pa) 

rc  characteristic  shear  stress  in  Ree-Eyring  model  (Pa) 
t h shear  stress  when  viscosity  equals  in  Ellis  model  (Pa) 
rm  shear  stress  at  conduit  center  (Pa) 

Tr  shear  stress  at  pipe  wall  (Pa) 
tw  shear  stress  at  conduit  wall  (Pa) 

spatial  domain  in  Dirichlet  functional  (m) 

B slit  half  thickness  (m) 

C viscosity  coefficient  in  Herschel-Bulkley  model  (Pa.sn) 

C viscosity  coefficient  in  Bingham  model  (Pa.s) 
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D Dirichlet  functional  (Pa2.m  x) 

/ Am7™ 

2 Fi  hypergeometric  function 

g 1 + / 

G constant  in  the  reduced  momentum  equation  (Pa.m-1) 

/p  Ca  definite  integral  for  Carreau  model  pipe  flow  (Pa3.s_1) 

IpCr  definite  integral  for  Cross  model  pipe  flow  (Pa3.s_1) 

JS)ca  definite  integral  for  Carreau  model  slit  flow  (Pa2.s_1) 

Is,Cr  definite  integral  for  Cross  model  slit  flow  (Pa2.s_1) 
k viscosity  coefficient  in  power  law  model  (Pa.s") 

K viscosity  coefficient  in  Casson  model  (Pa.s) 

L conduit  length  (m) 

m indicia!  parameter  in  Cross  model 

n flow  behavior  index  in  power  law,  Carreau  and  Herschel-Bulkley  models 
A p pressure  drop  across  conduit  length  (Pa) 

Q volumetric  flow  rate  (m3.s_1) 

r radius  (m) 

R pipe  radius  (m) 

s spatial  coordinate  representing  r for  pipe  and  z for  slit  (m) 

T stress  domain  in  Dirichlet  functional  (Pa) 
v fluid  speed  in  the  flow  direction  (m.s-1) 

W slit  width  (m) 

z spatial  coordinate  of  slit  thickness  (m) 

DM  variational  method  based  on  applying  Dirichlet  principle 
EL  variational  method  based  on  applying  Euler-Lagrange  principle 
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